Effect of long-range Coulomb interaction on shot-noise suppression 

in ballistic transport 
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We present a microscopic analysis of ^hot-noise suppression due to long-range Coulomb inter- 
action in semiconductor devices under ballistic transport conditions. An ensemble Monte Carlo 
simulator self-consistently coupled with a Poisson solver is used for the calculations. A wide range 
of injection-rate densities leading to different degrees of suppression is investigated. A sharp ten- 
dency of noise suppression at increasing injection densities is found to scale with a dimensionless 
Debye length related to the importance of space-charge effects in the structure. 



The phenomenon of shot noise, associated with the 
randomness in the flux of carriers crossing the active re- 
gion of a device, has become a fundamental issue in the 
study of electron transport through mesoscopic devices. 
In particular, the possibility of shot-noise suppression has 
recently attracted a-, lot of attention, both theoretically 
and experimentally.!!! At low frequency (small compared 
to the inverse transit time through the active region) the 
power spectral density of shot noise is given by Si — 
j2ql, where I is the dc current, q is the electron charge, 
and 7 is the suppression factor. When the carriers cross- 
ing the active region are uncorrelated, full shot noise with 
7=1 (Poisson statistics) is observed. However, correla- 
tions between carriers can reduce the shot-noise value, 
giving 7 < 1. In real mesoscopic devices different types 
of mechanisms resulting in shot-noise suppression can be 
distinguished: (i) statistical correlations due to the Pauli 
exclusion principle (important for degenerate materials 
obeying Fermi statistics) , (ii) short-range Coulomb inter- 
action (electron-electron scattering) , and (iii) long-range 
Coulomb interaction (by means of the self-consistent elec- 
tric potential). While the first two mechanisms Jiavc 
been extensively discussed in solid-,state literaturejil the 
last one has received less attention^ although its role in 
shot-noise suppression has been known for a long time in 
vacuum-tube devices.u The only exception that should be 
mentioned is the Coulomb blockade in resonant-tunneling 
devices, which can be also referred to the last mech- 
anism of suppression. The blockade is provided by a 
built-in charge inside a quantum well which redistributes 
the chemical potential and prevents the incoming car- 
riers from passing through the well, thereby resulting in 
carrier correlation and shot-noise suppression (see the ex- 
perimental evidencea). The Coulomb blockade is a con- 
sequence of long-range Coulomb interaction, and it acts 
under the sequential tunneling regime of carrier trans- 
port. 

The main objective of the present paper is to prove the 
importance of long-range Coulomb interaction between 



the carriers on the shot-noise power spectrum under the 
ballistic regime of electron transport. The ballistic regime 
is now accessible in modern mesoscopic devices like elec- 
tron waveguides, quantum point contacts, etc., which 
have characteristic lengths of the order, or smaller, than 
the carrier mean free path. The current existing theories 
invoked to interpret the experimentally observed shot- 
noise suppression in such devicesB~l3 assume that carriers 
move inside the device without inducing any redistribu- 
tion of the electric potential. We use a more rigorous 
approach which includes long-range Coulomb interaction 
between the carriers by considering the carrier transport 
in the self-consistent potential governed by the Poisson 
equation. We show that under the ballistic regime this 
interaction is crucial, and the noise characteristics are 
strongly modified depending on whether the carrier cor- 
relation (mediated by the field) is taken into account or 
not. 

To this purpose we consider a simple structure: a 
lightly doped active region of a semiconductor device 
sandwiched between two heavily doped contacts injecting 
the carriers into the active region. The device then acts 
similarly to a vacuum diode, with a relevant difference in 
the fact that there are two opposing currents instead of 
a single current. Electrons are emitted from the contacts 
according to a thermal-equilibrium Maxwell-Boltzmann 
distribution, and they move ballistically inside the active 
region (the mean free path is considered to be much larger 
than the distance L between the contacts) according to 
the semiclassical equations of motion. The fluctuating 
emission rate at the contacts is taken to follow a Poisson 
statistics. This means that the time between two con- 
secutive electron emissions is generated according to the 
probability density P(t) — re~ rt , where T = ^n c v t hS 
is the injection-rate density, with n c the electron density 
at the conta ct, S the cro ss sectional area of the device, 
and Vth = \/2ksT '/ (win) the thermal velocity (T is the 
lattice temperature, kg Boltzmann constant, and m the 
electron effective mass). The electron gas is assumed 
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to be nondegenerate to exclude possible correlations due 
to the Fermi statistics. For simplicity, the conduction 
band of the semiconductor is considered to be spheri- 
cally parabolic. Both the time-averaged current and the 
current fluctuations inside the active region of the device 
are analyzed for different bias voltages applied between 
the contacts. The calculations are performed by using an 
ensemble Monte Carlo simulator self-consistently coupled 
with a Poisson solver (PS). By using this approach we can 
analyze much more general situations than those studied 
in previous analytical calculations .□ 

For the calculations we used the following set of pa- 
rameters: T=300 K, m — 0.25mo, dielectric constant 
e = 11.7eo, sample length L = 2000A, and contact dop- 
ing n c ranging between 10 13 cm -3 and 4 x 10 17 cm -3 (al- 
ways at least two orders of magnitude higher than the 
sample doping). However, we must stress that the results 
we are going to present do not depend on the particular 
values of these parameters but only on the dimcnsionlcss 
length A = L/Ldc, where Ld c — \J eksT / q 2 n c is the De- 
bye length corresponding to the carrier concentration at 
the contact. 

Let us discuss briefly the steady-state spatial distribu- 
tions of the quantities of interest inside the sample. In a 
general case the carrier concentration is nonuniform, hav- 
ing maximum values at the contacts due to the electron 
injection and decaying toward the middle of the sample. 
Accordingly, without an external voltage bias the poten- 
tial distribution has a minimum in the middle of the sam- 
ple due to the space charge. When a positive voltage is 
applied to the anode, the minimum is displaced toward 
the cathode, while its amplitude tends to diminish. This 
minimum provides a potential barrier for the electrons 
moving between the contacts, so that a part of the elec- 
trons, not having enough energy to go over the barrier, 
are reflected back to the contacts. The most important 
fact is that the transmission through the barrier is cur- 
rent dependent, which is crucial in calculating the noise 
characteristics. In the structure the current is limited 
by the space charge and increases linearly with the ap- 
plied voltage up to a certain value of the external bias 
when the barrier vanishes so that all the electrons emit- 
ted from the cathode can reach the anode. Under the 
latter regime the cuirent is saturated and becomes inde- 
pendent of the bias.B It is important to stress that in our 
approach we do not impose a fixed number of electrons 
N to be present inside the sample. The value of N is 
determined by the emission rates of the contacts and the 
applied bias. Therefore, N fluctuates in time and we can 
evaluate both the time-averaged value (N) and its fluc- 
tuations by means of the Monte Carlo algorithm. One 
can observe that (TV) is constant at the increasing part of 
the current-voltage characteristic, and it decreases with 
the bias once the current is saturated. 

Under a fixed applied voltage the current density in 
the structure is given by I(t) — (q/L) J2fJi v i{t), where 
Vi(t) is the instantaneous velocity component along the 



field direction of the ith particle.Q The current auto- 
correlation function Cj{t) = (5I(t')8I(t' + t)) is evalu- 
ated from the sequence I(t) obtained from the Monte 
Carlo simulation, where the current fluctuation is given 
by SI(t) = I(t) - (I), To clarify the role of different 
contributions to the current noise we decompose the cur- 
rent autocorrelation function into three main contribu- 
tions Cj(i) = Cy(t) + CN(t) + CvN(t) associated, re- 
spectively, with the fluctuations in the mean velocity of 
electrons Cy, the fluctuations in the carrier number Cm, 
and the velocity- number cross correlation Cvn ■ The cor- 
responding formulas are given bycJ 

C v (t) = ^-(N) 2 (Sv(t')5v(t' + t)) (la) 

C N (t) = ^(v) 2 {6N(t')6N(t' + t)) (lb) 

Cvn® = £(v)(N){5v(t')6N(t' + t) 

+ SN(t')Sv(t' +*)> (lc) 

Fig. 1 shows the low-frequency value of the spectral 
density of current fluctuations Si = 2 Ci(t)dt nor- 
malized to 2gl,, where I s = qT = ^qn c VthS is the sat- 
uration current (notice that this is the maximum cur- 
rent that a contact may provide). This normalization is 
performed in order to compare the results for different 
injection-rate densities (different contact dopings). We 
provide the results for two different simulation schemes. 
The first one involves a dynamic PS, which means that 
any fluctuation of space charge appeared due to the ran- 
dom injection from the contacts causes a redistribution 
of the potential, which is self-consistently updated by 
solving Poisson equation at each time step during the 
simulation to account for the fluctuations associated to 
long-range Coulomb interaction. In the second scheme 
we use a static PS to calculate only the stationary poten- 
tial profile, and once the steady state is reached, the PS is 
switched off, so that the carriers move in the frozen non- 
fluctuating electric field profile. We have checked that 
both schemes give exactly the same steady-state spatial 
distributions and average current, but the noise charac- 
teristics are different. Several values of n c (and therefore 
several injection-rate densities) have been considered. As 
n c increases, space-charge effects become more and more 
significant, the dimensionless parameter A being the in- 
dicator of their importance. 

In the static case, by increasing the applied voltage 
U we always obtain, an excellent coincidence with the 
well-known formulatLj used to describe the crossover from 
thermal to shot noise when carrier correlation play no role 
(represented in the figure by dashed lines): 

Si = 2q = 2qIcoth{qU/2k B T) (2) 

where I = I + — I~ is the total current flowing through 
the diode, consisting of two opposing currents, I + = 
/ s exp[— qV m /kBT] in the forward-bias direction and 



2 



I~ = I s exp[— q(V m + U)/IcbT] in the opposite direction, 
V m being the potential minimum induced by the space 
charge, which is dependent on U . This agreement sup- 
ports the validity of the simulation scheme used for the 
calculations. For qU <C fcgT, I + ~ I~, thermal noise 
is dominant and Si ~ AqI s cxp[—qV m /kBT]. Therefore, 
for the lowest value of A (when space charge is negligi- 
ble and V m — > 0), Si — > iql s , while as A increases V m 
becomes significant and Si decreases. When qU > fc^T, 
I + 3> I~ ■, the transition from thermal noise to shot noise 
takes place and Si ~ 2ql + . Finally, for the highest values 
of U, saturation occurs, V m vanishes, and Si » 2ql s . 

For the lowest values of A no difference between the dy- 
namic and static cases is obviously detected. However, 
for higher A, when space-charge effects become signifi- 
cant, the picture is drastically different for the dynamic 
case. Starting from qU ~ ksT the current noise, instead 
of increasing, decreases until the proximity of saturation. 
Under saturation, the results for both schemes coincide 
and full shot noise 2ql s is recovered. When compared 
with the static case the noise suppression is stronger for 
higher A (more important space-charge effects) . 

To understand the physical reason for the shot-noise 
suppression, in Fig. 2 we provide the decomposition of Si, 
calculated with static and dynamic PS, into the additive 
contributions Sv, Sn, and Svn [Eqs. (1)] for different ap- 
plied voltages U and A=7.72 (n c = 2.5 x 10 16 c m - 3 ). The 
contributions of Sn and Svn to the current noise vanish 
at equilibrium (U — > 0), since they are proportional to 
(v) 2 and (v) — > 0. Thus, for small biases (qU <C k B T) 
Si ~ Sy, which means that the current noise is thermal 
noise associated with velocity fluctuations and is gov- 
erned by the Nyquist theorem Si « 4ksTG, with G — 
dI/dV\v=o the conductance. For this case the results for 
the static [Fig. 2(a)] and dynamic [Fig. 2(b)] schemes ev- 
idently coincide. However, starting from qU ~ fc^T the 
difference becomes drastic. For the dynamic case the 
velocity- number correlations, represented by Svn, are 
negative, while for the static case they are positive. Fur- 
thermore, for the current fluctuations calculated using 
the self-consistent potential, Sn and Svn are of opposite 
sign and compensate for each other, so that Si approxi- 
mately follows Sv as long as the current is space-charge 
limited. As a consequence, the current noise, which now 
corresponds to shot noise, is considerably suppressed be- 
low the value 2ql given by the static case. This result 
reflects the fact that as the carriers move through the ac- 
tive region, the dynamic fluctuations of the electric field 
modulate the transmission through the potential mini- 
mum and smooth the current fluctuations imposed by the 
random injection at the contacts. Therefore, the coupling 
between number and velocity fluctuations induced by the 
self-consistent potential fluctuations is the main respon- 
sible, through Svn, for the shot-noise suppression. This 
velocity-number coupling becomes especially pronounced 
just before the current saturation (U ks 7ksT/q), when 
the potential minimum is close to vanishing completely 
(V m — > 0) and the fluctuations of the potential bar- 



rier modulate the transmission of the more populated 
states of the injected carriers (the low- velocity states). 
Under saturation conditions space-charge effects do not 
modulate the random injection (no potential minimum is 
present) and again both dynamic and static cases provide 
the same additive contributions and total noise (2ql s ). 

Finally, in Fig. 3 we present the reduction factor 7' 
defined as the ratio between Sj as calculated with the 
dynamic PS (Si-d) and as given by Eq. @ neglecting 
the influence of long-range Coulomb interaction, 7' = 
Si-d/ftql coth(qU/2k B T)]. In the context of our calcu- 
lations 7' is more appropriate than the standard sup- 
pression factor 7 = Si-d/{2ql), since it covers both the 
thermal and shot-noise range of applied voltages. Here 
it is observed how the shot-noise suppression becomes 
more pronounced as A increases. For example, for A=30.9 
it reaches 0.04. Thus, our self-consistent approach pre- 
dicts much lower values of the suppression factor than the. 
previous analytical model of van der Ziel and Bosman,EI 
where the dependence of the potential minimum and its 
position on the applied voltage was not taken into ac- 
count. 

In principle, the value of the parameter 7' in our model 
has no lower limit. We observe that it follows asymp- 
totically the behavior 7' ~ ksT/qU in the range where 
shot-noise suppression is more pronounced (qU 3> fc^T, 
U <C U sa t). 7' can reach a value as low as desired by 
appropriate increasing the sample length and/or the car- 
rier concentration at the contact, provided the transport 
remains ballistic. However, with increasing the device 
length (or lattice temperature) the carrier transport ac- 
tually goes from the ballistic to the diffusive regime, and 
the shot-noise suppression is washed out. Therefore, the 
maximum suppression factor predicted by our calcula- 
tions for a system with a given value of the mean free path 
A p would be obtained approximately at A max ~ \ p /Ld c . 
Moreover, when the carrier concentration at the contact 
is increased so that the electron gas becomes degenerate, 
statistical (Pauli) correlations between the carriers ap- 
pear, which will be additive (in the sense of shot-noise 
suppression) to the Coulomb correlations. 

It should be emphasized that two essential conditions 
are necessary for the strong shot-noise suppression due 
to long-range Coulomb interaction: (i) the presence of a 
potential barrier inside the device which controls the cur- 
rent, and (ii) the carrier transmission through the barrier 
should depend on the current. This fact is quite gen- 
eral and, therefore, the results obtained in the present 
paper extend to much more physical situations. For ex- 
ample, in the recent experiment by Reznikov et aln the 
shot-noise level measured in a quantum point contact in 
the pinched-off regime was found to be unexpectedly low 
(about one-third and less). In that regime the trans- 
port is controlled by the potential barrier present at the 
gates, and the both conditions for the shot-noise sup- 
pression mentioned above are fulfilled. Hence the results 
obtained in our calculations strongly support the sugges- 
tion of the authors of the experiment that the origin of 
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the discrepancies between experimental results and theo- 
retical predictions lies in the disregarding of Coulomb in- 
teraction between electrons passing through the contact. 
More precisely, the electron flow considerably modifies 
the potential distribution inside the contact, yielding the 
coupling of velocity- number fluctuations, and resulting 
in shot-noise suppression. 

In conclusion, we have investigated the influence of 
long-range Coulomb interaction on shot-noise suppres- 
sion in ballistic transport by using an ensemble Monte 
Carlo simulator self-consistently coupled with a Poisson 
solver. We have found that this suppression is stronger 
as space-charge effects become more important, and it 
can be monitored by a dimcnsionless parameter A. More 
than one order of magnitude of shot-noise suppression is 
predicted. The main contribution to the suppression is 
found to originate from the velocity-number correlations 
induced by the self-consistent field. 

This work has been partially supported by the 
Comision Interministcrial de Ciencia y Tecnologfa 
through the project TIC95-0652. 



FIG. 1. Current-noise spectral density Si vs applied 
voltage U calculated by using static (open symbols) and 
self-consistent (close symbols, solid line) potentials for sev- 
eral injection-rate densities n c (in cm -3 , with the corre- 
sponding A): (o) 10 13 , A=0.15; (O) 2xl0 15 , A=2.18; (y) 
10 16 , A=4.88; (A) 2.5xl0 16 , A=7.72; (□) 10 17 , A=15.45; (0) 
4xl0 17 , A=30.9. The static case is shown to be nicely de- 
scribed by Eq. (2) (dashed line). The dotted lines represent 
2ql (marked for each injection-rate density by the correspond- 
ing symbol). 

FIG. 2. Decomposition of the spectral density of 
current fluctuations Si into velocity, number and veloc- 
ity-number contributions vs applied voltage for the case 
n c =2.5x 10 16 cm -3 , A=7.72 calculated by using (a) the static 
and (b) the dynamic Poisson solver. 

FIG. 3. Shot noise reduction factor 7' vs voltage U for 
several injection-rate densities n c (different values of A). 
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